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Q ■ Abstract 



\ For a 3-qubit Heisenberg model in a uniform magnetic field, the pairwise thermal 

' entanglement of any two sites is identical due to the exchange symmetry of sites. In this 

paper we consider the effect of a non-uniform magnetic field on the Heisenberg model, 
1^ \ modeling a magnetic impurity on one site. Since pairwise entanglement is calculated 

I I by tracing out one of the three sites, the entanglement clearly depends on which site 

^ ■ the impurity is located. When the impurity is located on the site which is traced out, 

^ , that is, when it acts as an external field of the pair, the entanglement can be enhanced 

qh| to the maximal value 1; while when the field acts on a site of the pair the corresponding 

concurrence can only be increased from 1/3 to 2/3. 

• I— I 

X. 

1 Introduction 

There is currently an ongoing effort to study entanglement in multipartite systems, since 
such entangled states may provide a valuable resource in quantum information processing 
Ijl]]. Recently entanglement in quantum operations p|, |] and in indistinguishable fermionic 
and bosonic systems ^ 0] have been considered. Entanglement in two-qubit states has 
been well studied in the literature. Various kinds of three-qubit entangled states have also 
been studied [l^ , which have been shown to possess advantages over two-qubit states 
in quantum teleportation |JlT[, dense coding [|1^] and quantum cloning |T3[] . 



One interesting and natural type of entanglement, thermal entanglement, was introduced 



and studied in the context of the Heisenberg XXX |T^, XX |T^, and XX Z [0] models as 



well as the Ising model in a magnetic field The Heisenberg interaction has been used 



to simulate a quantum computer [113, and can also be realized in quantum dots |TH], nuclear 
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electronic spins and optical lattices [^. By suitable coding, the Heisenberg 
interaction can be used for quantum computation |2^. Entanglement in the ground state of 
the Heisenberg model has been discussed previously In an earlier note [2|] we presented 
an analytical study of pairwise entanglement in the 3-qubit Heisenberg model in a uniform 
magnetic field and found that the magnetic field can greatly enhance pairwise entanglement. 
Due to exchange symmetry in this cyclic model the entanglement of any two sites is identical. 

In this paper we consider the effect of a magnetic impurity on entanglement in the 
Heisenberg model. We find unsurprisingly that the effect of such an inhomogeneous magnetic 
field on the entanglement depends on which site the impurity is located, although in a non- 
intuitive way. When the field may be considered as an external field of the pair, that is, 
when it is located on the site which is traced over, then it can enhance the entanglement to 
its maximal value, as measured by the concurrence. When the field acts on a site of the pair 
the concurrence can be increased from 1/3 to 2/3, but not to its maximal value 1. 



2 XX Heisenberg model with magnetic impurity 

We consider the 3-qubit XX Heisenberg model in a magnetic field acting on the third site 
only. The Hamiltonian is |2^ 



H 



T 3 



3' 



(1) 



1=1 



where we use BJ rather than B to denote the magnetic field. The Hamiltonian (p has eight 
distinct eigenvalues when i? 7^ 



-J5, 



E2 = - J(l + B) 
J 



E^ = ^il + B_] 



B 



+) 5 



-J(l-i?), 



B_ 



(2) 



^ 2 
B+) , E7 = JB, 

where B± = (45^ ±AB + 9^/^. When B = 0, the energy levels are degenerate 

El = Ej = 0, El = £^3 = i?4 = 2J, E2 = E^ = E^ = — J. (3) 

In the antiferromagnetic case (J > 0), the ground state is i?2, while in the ferromagnetic 
case (J < 0), the ground state is i?4. 

The corresponding non-degenerate, orthogonal eigenstates are 

00) = |000), 

01) =A/'i(|100) + |010) + ai|001)), 

02) = 2-1/2 (1010) -1100)), 

03) = 2-1/2 (1101) -1011)), 

04) =Ar4(a4|110) + |101) + |011)), 

05) =Ar5(|100) + |010) + a5|001)), 

06) =Ar6(a6|llO) + |lOl) + |Oll)), 

h) = |iii), 



(4) 
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where 



ai 
04 



2 2^ 



B, 



as 



2 2 ^ 



(5) 



and = (2 + a2)-i/2 = 1^ 4^ 5^ 5) are normalization constants. 
It is interesting to note that the eigenvalues transform under B 



-Bhy 



Eo^Ej, Ei^Ei, E2^Es, Er, ^ Eq, (6) 

and so the a^'s transform by ai ^ a^, ^ a^. This leads to invariance of the entanglement 
under B ^ —B. 

The density operator p(T) at temperature T can be written as 



Z 



i=0 



where j3 = 1/kT and Z is the partition function 
Z = tr(e-^^) =f]e-^^' 

= 2(1 + e^^) cosh{J(3B) + 26"^'^/' 



cosh -JPB+ + cosh -JPB 



1 



(7) 

(8) 
(9) 



3 Concurrence of pairwise entanglement 



The easiest way to calculate the entanglement is by means of the concurrence C [|26[ between 
a pair of qubits, which is defined as 



C = max {Ai — A2 — A3 — A4, 0} , 
where the quantities Aj are the square roots of the eigenvalues of the operator 

^ = p(or^®af)p*K®af) 



(10) 



in descending order; p is the density operator of the pair and it can be either pure or mixed. 
The entanglement of formation is a monotonic function of the concurrence C, varying between 
a minimum of zero for C = 0, and a maximum of 1 for C = 1. 

We now derive the concurrence for any pair of sites in our model. Due to symmetry 
under the exchange of sites 1 and 2, the entanglement between sites 1 and 3 is the same as 
that between sites 2 and 3, and so we need only consider entanglement between sites 1 and 
3, and between sites 1 and 2. 

Taking the trace over the second (third) site, we can obtain the reduced density operator 
Pi3(pi2) of the sites 1 and 3 (1 and 2). Both and pis take the following form 



/ u 



1 
Z 



Wi 

y 



y 

W2 



\ 



V J 



(12) 
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Here, for pi2, the nonzero matrix elements are given by 



e 



V 



2 2 



2 2 

e-^^o + a^AT^e-^^i + a2Ar2e-^^% 

e-'^^^ + a^AT^e-^^^ + a^'e-^''^ (13) 



while for the pi3 case, we have 

w, = alMle-"''^ + ^e-'^^^ + AT^e"^^^ + a^A^^e"^^^ + AT^e'^^^ 

e-^^o + ATfe-^^^ + ^e-'^^^ + A^2e-^^^ 
1 



e 



+ - -/3i?3 + ^Ule-"^^ + AT^e-^^^ (14) 



2 

The concurrence has the form 

2 r ^ 

C = — max ||?/| — vtiw, 0| . (15) 

The system is entangled when C > 0, and maximally entangled when C = 1. The exchange 
interaction constant J and the temperature T always appear in the form J/ kT in the con- 
currence and thus we can define the scaled temperature t = kT/\J\ > 0. The concurrence is 
a function of r and B. 

From Eqs. (§) it is easy to see that y ^ y and u v when B — ^ —B. This means that 
the concurrence is invariant under B <-> —B; 

C{t,B)=C{t,-B). (16) 

We therefore only consider the case B > case hereafter. 



4 Discussion and results 
4.1 Ci2 

We first consider the entanglement between sites 1 and 2. In Fig. 1 and 2 we give plots of 
the concurrence of pu against r and B. We know that entanglement appears only in the 
antiferromagnetic case (0 < r < 1.27) when 5 = p4[ (also see Figure 1). From Fig. 
1 and 2 we see that, when the magnetic impurity is located on the third site, both the 
antiferromagnetic and ferromagnetic cases are entangled in the range < r < tq, where tq 
depends on B. 
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Fig. 1 also suggests that the concurrence Cu goes to 1, namely, the sites 1 and 2 reach 
maximal entanglement, when r ^ for large enough B, in both the antiferromagnetic and 
ferromagnetic cases. This fact can be shown analytically as follows. 

Consider first the antiferromagnetic case (J > 0). In this case E2 is the ground state; 
that is, E2- Ei<0 for alH 7^ 2 and thus e"^^^ ^ ^-PE, ioi i ^ 2 in the hmit r ^ 0. Note 
that all A/i and are finite. Then we have 

2/-^e-^^^ Z^e-^^\ 1-0, 1-0, (17) 

namely, C12 — 1 when r — 0. 

For the ferromagnetic case (J < 0), one can check that E4 — Ei < for all i ^ A and 
g-/3E4 ^ ^-PEi ^ - _^ ^Yie hmit r — 0. Then we have 

y^Af!e-^^\ Z-e-^^^ | - 0, | - a^AT^, (18) 
namely, 

C12 - 2Ar,2 = (19) 

when r — s> 0. In the limit 5 — > 00, 04 — > and therefore C12 — 1. In the limit 5 — +0, but 
5 > r, C12 — 2/3. 



It is interesting to note that, when B = 0, C12 — 1/3 in the limit r — In this 



case the ground state is 3-fold degenerate and the approximation we used above is not valid. 
This again indicates the role of degeneracy in entanglemant. 



4.2 C 
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We consider the entanglement between sites 1 and 3. /^From Fig. 3 and 4 we see that: 

1. In contrast to the 1-2 case, the concurrence increases to a maximum with increasing 
B and then decreases. The lower r, the smaller B at which the concurrence reaches its 
maximum value. 

2. For small B, entanglement occurs only in the ferromagnetic case (J < 0), while 
for large enough B (e.g. B = 10), entanglement occurs in both the antiferromagnetic and 
ferromagnetic cases, but it is very weak. 

Fig. 4 suggests that the maximal entanglement occurs in the ferromagnetic case when 
r — and B is also much smaller than 1. In this case, Ei is very close to the ground state 
E4 and exp^—PE^) and exp{—i3Ei) are much bigger than others. We can also check that 

exp ( > 1, (20) 



and that 



Afi ~ Ma 



1 

3' 



Ctl ~ ^4 ~ 1. 



The concurrence is then given approximately by 



C 



13 



exp 


(ID 




1 + exp 


(If) 





(21) 



(22) 
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from which we conclude that the maximal concurrence is 2/3 when B is much bigger than 
r and much smaller than 1. 



In summary, we list our results in the following table. 





Maximal concurrence 


Entanglement ranges 


B=0 


1/3 


Antiferromagnetic case only 


12 


1, when r ^ and 
B is big enough. 


In both ferromagnetic and antiferromagnetic 

cases 


13 


2/3 

for antiferromagnetic 
case and r <^ B ^ 1 


When B is small, only the ferromagnetic case is 
entangled. When B is big enough, both 
the antiferromagnetic and ferromagnetic cases 
are entangled, but the entanglement is very weak. 



5 Conclusion 

In this paper we considered the effect of a non-uniform magnetic field on the Heisenberg XX 
model, modeling a magnetic impurity on only one site. In contrast to the uniform magnetic 
field case p4[ where the pairwise thermal entanglement of any two sites is identical due 



to the exchange symmetry of sites, the entanglement due to a non-uniform magnetic field 
clearly depends on which site the impurity is located. When the impurity is located on the 
site which is traced out, that is, when it acts as an external field of the pair, the concurrence 
corresponding to the entanglement can be enhanced to the maximal value 1 from 1/3; while 
when the field acts on a site of the pair the concurrence can only be increased from 1/3 to 
2/3. Maximal entanglement is achieved when the temperature tends to zero. 



In [£4[, the entanglement was related to the degeneracy of the system. In the present 
model, the magnetic field removes all the degeneracy of the energy levels present when B = 
and the entanglement is thus greatly enhanced. 
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Figure 1: Concurrence C12 against r for different magnetic field B — 0,1,10. 
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Figure 2: Concurrence C12 against B at different temperature r = 0.1, 0.5 and 1. 
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Solid line: Ferromagnetic case 
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Figure 3: Concurrence C13 against r for different magnetic field B = 0,1, 10. For antifer- 
romagnetic case (dotted line) with B — 10, the entanglement occurs although it is very 
week. 
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Figure 4: Concurrence C13 against B at different temperature. For antiferromagnetic case 
(dotted line), r — 2. 
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